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ABSTRACT 

Isolated low-mass stars are formed, in the standard picture, from the collapse 
of dense cores condensed out of strongly magnetized molecular clouds. The dy- 
namically collapsing inflow traps nearly half of the critical magnetic flux needed 
for the core support and deposits it in a small region surrounding the proto- 
star. It has been argued previously that the deposited flux can slow down the 
inflow, allowing matter to pile up and settle along field lines into a magnetically 
supported, circumstellar disk. Here we show that the disk typically contains 
~ 10% of the stellar mass, and that it could become self-gravitating under plau- 
sible conditions during the rapidly accreting, "Class 0" phase of star formation. 
Subsequent fragmentation of the self-gravitating, magnetically subcritical disk, 
driven by magnetic diffusion, could produce fragments of substellar masses, which 
collapse to form brown dwarfs and possibly massive planets. This scenario pre- 
dicts substellar object formation at distances of order 100 AU from the central 
star, although orbital evolution is possible after formation. It may provide an ex- 
planation for the small, but growing, number of brown dwarf companions found 
around nearby stars by direct imaging. The relatively large formation distances 
make the substellar companions vulnerable to dynamic ejection, particularly in 
binary (multiple) systems and dense clusters. Those ejected may account for, at 
least in part, the isolated brown dwarfs and perhaps free-floating planetary mass 
objects. 

Subject headings: accretion, accretion disks — ISM: clouds — MHD — stars: 
formation — stars: low-mass, brown dwarfs 

1. Introduction 

The past several years have witnessed remarkable progress in the study of brown dwarfs 
(Basri 2000). To date, more than 200 of such objects have been discovered. Most are found 
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in isolation; only a small fraction orbit around stars, typically at separations of tens of AUs 
and beyond. Despite the impressive pace of discovery, investigation into their origins has only 
just begun. The isolated brown dwarfs could be formed in small multiple systems through 
cloud fragmentation like stars, but ejected through interaction with other (heavier) members 
before gaining enough mass for hydrogen burning (Reipurth & Clarke 2001). The brown 
dwarf companions could in principle be produced in (rotationally supported) circumstellar 
disks through gravitational fragmentation, although it is unclear whether the conditions 
for fragmentation can arise naturally during star formation (Durisen 2001). In this Letter, 
we advance a new scenario for forming substellar companions, through the gravitational 
fragmentation of magnetically, rather than rotationally, supported disks. Such disks are 
expected to develop around protostars formed out of strongly magnetized clouds (Li & 
McKee 1996; LM96 hereafter), as envisioned in the standard picture of isolated low-mass 
star formation (Shu, Adams & Lizano 1987). 

Conceptually, the formation of magnetically supported disks parallels that of the more 
familiar rotationally supported disks. Once gravity has initiated the dynamic collapse of a 
rotating magnetized dense core of molecular cloud, both the magnetic flux and the angular 
momentum associated with rotation are trapped by the collapsing flow and carried to the 
vicinity of the protostar. However both must eventually be stripped almost completely 
from the matter before it enters the star; otherwise, the stellar magnetic field strength and 
rotation rate would be many orders of magnitude higher than actually observed. These 
are, respectively, the magnetic flux and angular momentum problem of star formation. The 
stripping of angular momentum is thought to take place primarily in a rotationally supported 
disk by friction between adjacent annuli of matter moving at different speeds. The mechanism 
for magnetic flux stripping is less clear. LM96 proposed that it occurs as the field lines 
decoupled from the stellar matter expand against the collapsing inflow, reversing the inward 
motion of charged particles and the magnetic field tied to them. The bulk neutral material 
can still slip through, but only at a much reduced speed because of frequent collisions with 
the already stopped charged particles. The slowdown causes the infalling matter to pile up 
and settle along field lines into a flattened structure — a magnetically supported circumstellar 
disk (see also Ciolek & Konigl 1998 and Contopoulos, Ciolek & Konigl 1998). 

We show in this Letter that the magnetically supported disk could become self-gravitating 
under plausible conditions during the earliest, "Class 0" phase of star formation, when the 
accretion rate is highest (§ 2). The self-gravitating disk is magnetically subcritical, and 
fragmentation driven by magnetic diffusion could lead to substellar object formation (§ 3). 
We discuss the expected characteristics of the substellar companions and their observational 
implications in § 4. 
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2. Self-Gravitating Magnetically Supported Disks 

We consider the simplest case of axisymmetric non-self-gravitating disk and derive a 
condition for the disk to become self-gravitating around a single, isolated star M*; possible 
nonaxisymmetric effects will be commented upon toward the end of the section. We estimate 
the disk properties in a way similar to LM96, but include the magnetic compression, which 
controls the disk density, and the magnetic tension force, which dominates the magnetic 
pressure gradient in the radial force balance: 

GMjr 2 « B z B r /2nE. (1) 

We show below that the disk column density S is roughly independent of the cylindrical 
radius r. To counter the larger stellar gravity at a smaller radius, the field strength must 
increase inward, roughly as r _1 . For such a distribution, the radial and vertical field com- 
ponents are comparable on the disk (assuming a potential field outside as usual). They 
are 

B z m B r Ki (27rGSM,) 1/2 /r. (2) 
The magnetic flux enclosed within a disk of radius R is 

r R 

$ d = / 2nB z rdr « 2irR(2irGEM^) 1/2 . (3) 
Jo 

Since this flux is mostly that stripped from the stellar mass M* (LM96), we have $d « = 
€(2ttG 1 ^ 2 )M^, where e (< 1) is the flux-to-mass ratio of the dynamically collapsing envelope 
that feeds the star-disk system in units of the critical value 2-nG 1 / 2 . Eliminating $ d from 
equation (3), we obtain the following estimate of the disk mass 

M d « 7ri? 2 S « 0.5e 2 M*, (4) 

which is simply the mass needed to "weigh down" the magnetic flux released from the star 
and prevent it from escaping. Even though the star-disk system remains magnetically su- 
percritical, the disk (which contains a small fraction of the mass of the system but most of 
its flux) is subcritical, by a factor of ~ 2/e. Calculations of the evolution of strongly mag- 
netized clouds driven by ambipolar diffusion suggest that e is close to half (e.g., Nakamura 
& Li 2002). This value would yield a disk mass of ~ 10% of the stellar mass, or some 10 2 
Jupiter masses for solar-mass stars. Therefore, there appears to be enough matter in the 
disk to form one brown dwarf near the hydrogen burning limit, a few less massive brown 
dwarfs, or several massive planets. 

The mass density p of the disk is determined by balancing the outward thermal pressure 
gradient in the vertical direction against the squeezing of the stellar tidal force, disk self- 
gravity, and the radial component of the magnetic field. It is easy to show that the magnetic 
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contribution dominates and 

p « B 2 /8na 2 « GM*£/4a 2 r 2 , (5) 

where a = 0.188 T^ 2 km s" 1 is the isothermal sound speed and Tio = T/10 K. Equation (2) 
was used to derive the second relation. The column density X at any radius r is related to 
the local mass accretion rate M and infall speed V through £ = M/2nrV. Approximately, 
V is the speed with which the bulk neutral matter slips across the field lines, which remain 
more or less fixed in space (LM96). The slippage speed can be written as a product of the 
total magnetic force per unit mass (ss GM*/r 2 ) and the magnetic coupling time t c . The 
latter is related to the local free-fall time tg = (3n/32Gp) 1 ^ 2 through t c = tg/ts, where 
the characteristic flux loss time ts from a magnetically supported cloud has been discussed 
extensively by Nakano and collaborators. According to Nakano, Nishi & Umebayashi (2002), 
ts is a factor of ~ 10 — 10 2 larger than tg up to a density of ~ 10 10 — 10 11 cm -3 for a cosmic ray 
ionization rate of ~ 10~ 17 — 10~ 16 s _1 and an MRN grain size distribution (see their Fig. 3). 
It drops below ts beyond ~ 10 12 cm -3 , as both ions and charged grains become decoupled 
from the field lines. The decoupling density depends somewhat on the grain size distribution. 
It could be increased significantly due to rapid grain growth and/or sedimentation toward 
the disk midplane (Sano et al. 2000). We will concentrate on the part of the disk that is 
well coupled to the magnetic field (which makes mass pileup possible), and adopt t B = 20 t s 
for a rough estimate. In this case, t c = iff/20 and the column density is given by 

S « 9 M 2 /a 2 M, = 5 rtfm^T- 1 g cm" 2 , (6) 

where m is the accretion rate normalized by 10~ 5 M Q yr _1 and m* the stellar mass in M & . As 
advertised, the column density is roughly independent of radius r as long as the accretion rate 
does not vary much over the disk and the disk is more or less isothermal. For a canonical 
dust opacity of 0.01 cm 2 g^ 1 , the disk is typically optically thin, making the isothermal 
approximation justifiable. 

In the absence of a strong magnetic field, a circumstellar disk becomes gravitationally 
unstable when the radius of a Jeans-mass clump is smaller than the Roche value. This 
condition roughly corresponds to the Toomre's criterion (e.g., Shlosman & Begelman 1989) 

Q = afi/vrGS £ 1, (7) 

where f2 = (GM^/r 3 ) 1 / 2 is the Keplerian frequency, even though the disk is not rotationally 
supported. The presence of a strong magnetic field does not suppress the Jeans instability 
in the lightly ionized disk; it merely lengthens the growth timescale (Langer 1978; see § 3 
for further discussion). Near the outer edge of the disk 

R w (M d /7rS) 1/2 w 2.6 x lOh^m^m^ 2 AU, (8) 



- 5 - 



where most of the disk mass is located, the Toomre Q parameter has the approximate value 

Q R « 0.8 e -f m-^Tif , (9) 

where e 5 is the flux-to-mass ratio e normalized by the typical value 1/2. The most uncertain 
quantity in the above expression is the accretion rate m. If it is large enough (m ^ 1), the 
disk could become self-gravitating (i.e., Qr < 1). 

Low-mass protostars are thought to grow most rapidly during the so-called "Class 0" 
phase (Andre, Ward-Thompson & Barsony 2000), which lasts for only a few times 10 4 years. 
To accrete the bulk of the mass of Sun-like stars over such a short period of time, the accretion 
rate must be of order 10~ 5 M yr _1 or higher. Such high accretion rates are inferred from the 
observations of molecular outflows from Class sources (Bontemps et al. 1996), and found in 
core collapse calculations (e.g., Tomisaka 1996; Li 1998). We conclude that during the Class 
phase of star formation the magnetically supported disk could become self-gravitating and 
thus prone to fragmentation into small pieces. Before discussing fragmentation, we comment 
on a possible complication. 

The disk has a flux-to-mass ratio B z /H increasing toward the center, and should be 
unstable to (nonaxisymmetric) interchange instability (LM96; Ciolek et al. 1998) according 
to the criterion of Spruit & Taam (1990), derived in the limit of frozen-in magnetic field 
and no self-gravity. Nonlinear developments of the instability could potentially drain the 
disk material faster than estimated above, and thus prevent the disk from becoming self- 
gravitating. However, in our case the magnetic field is not frozen in the matter and the disks 
relevant to companion formation are self-gravitating. Moreover, a realistic disk would have 
a significant amount of differential rotation and be magnetically linked to the surrounding 
medium, both of which tend to stabilize the instability (Spruit, Stehle & Papaloizou 1995), 
as does the (likely differential) twisting of the disk field lines associated with magnetic 
braking (Sakurai 1989). Nonaxisymmetric calculations are needed to ascertain whether the 
combination of magnetic diffusion, self-gravity, rotation, and field line linkage and twisting 
can suppress the instability altogether and, if not, to what extent nonlinear developments of 
the instability modify the disk structure. 



3. Magnetic Diffusion-Driven Disk Fragmentation 

The disk is magnetically subcritical. If frozen in matter, the magnetic field would 
prevent fragmentation completely (Nakano 1988). However, in a lightly ionized medium, the 
presence of a strong magnetic field does not change the criterion for fragmentation; rather, 
the fragmentation occurs on the (longer) magnetic diffusion, instead of dynamic, timescale 
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(Langer 1978). The minimum mass of fragments will still be given by the (thermal) Jeans 
mass, which in a flattened geometry has the following approximate expression (Larson 1985) 

Mj = 1.17a 4 /G 2 H Ri3x lO~ 3 m~ 2 m*T? M . (10) 

It corresponds to ~ 3 times the Jupiter mass (M Jup ) for typical parameters, which is in the 
planetary mass range. This typical Jeans mass is smaller than the oft-quoted opacity-limited 
minimum value (~ 7 Mj up ; Low & Lynden-Bell 1976; Silk 1977), although in a disk geometry 
where radiation can escape more easily this minimum could be lower. 

Without magnetic field, a self-gravitating fragment of one Jeans mass will contract in a 
shorter time onto itself than toward the central star. This remains true even in the presence 
of a strong magnetic field, since both timescales are lengthened by the same factor, set by 
magnetic diffusion. The timescale for the magnetic field to diffuse out of a fragment of 
Jeans radius Rj = (Mj/ttS) 1 / 2 is set roughly by r R3 Rj/V, where V « t c (GMj/R 2 ) is 
the slippage speed of neutral matter across field lines driven by self-gravity. Adopting the 
same approximation for the magnetic coupling time t c = tg/20 as before, and computing the 
freefall time t$ near the disk outer edge where 

P H ~3x 10- 14 e- 2 m 4 m; 2 T- 3 g cm" 3 (11) 

[obtained from equations (5), (6) and (8)], we find 

r^3x 10 4 e7 5 1 m- 2 m*T 3 /2 yrs, (12) 

which for typical parameters is comparable to the lifetime of Class phase. Once the 
fragment has become magnetically supercritical, it collapses dynamically from inside out to 
form one or more substellar objects, as in the standard scenario for forming isolated low-mass 
stars through ambipolar diffusion (Shu et al. 1987). The main difference lies in the density 
of the background, magnetically subcritical material, which is of order ~ 10~ 20 g cm -3 in the 
stellar case but ~ 10~ 14 g cm -3 in the substellar case. The density, of course, sets the mass 
scale (for a ~ 10 K gas). In the substellar case, the density is high enough that magnetic 
decoupling could be a potential problem; its complete treatment is hindered by uncertainties 
in the grain properties of the disk. 

We therefore arrive at the following plausible scenario for the production of substellar 
companions, based on the standard picture of isolated low-mass star formation. To avoid 
the famous "magnetic flux problem", most of the magnetic flux originally trapped by the 
collapsing core material that enters the star must be stripped and deposited in a magneti- 
cally supported disk containing ~ 10% of the stellar mass. If the disk material is not drained 
at too large a speed by interchange instability and/or magnetic decoupling, it could achieve 
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a high enough surface density during the Class phase to become self-gravitating. Subse- 
quent contraction of self-gravitating pockets of substellar masses across field lines, over a 
magnetic diffusion timescale comparable to the Class lifetime, could produce magnetically 
supercritical fragments, which collapse to form substellar companions. This scenario does 
not involve rotation directly, and should not be affected by density waves, which tend to 
suppress gravitational fragmentation of rotationally-supported disks (Durisen 2001). 

4. Implications for Brown Dwarf Companions 

If the substellar companions are formed from the fragmentation of magnetically sup- 
ported disks, they should have initial masses between the opacity-limited Jeans mass, which 
is in the (upper) planetary mass range, and the disk mass, typically comparable to that of 
a massive brown dwarf. Further mass increase by accretion from the residual envelope is 
possible, but probably not by much because of rapid decline of the overall accretion rate 
after the Class phase, as a result of, e.g., the powerful protostellar wind unbinding much of 
the envelope. If substantial mass accretion does occur, the substellar objects may be turned 
into (very) low mass stellar companions. 

A prediction of our scenario is that the substellar companions should form at distances 
of a few hundred AUs from the central star, set by the typical size of the self-gravitating disk 
(cf. equation 8). Interestingly, more than half of the brown dwarf companions detected so far 
have separations in the range of ~ 50 — 300 AU (see Reid et al. 2001 for a compilation). The 
orbits could shrink after formation for several reasons: first, the disk material out of which 
the objects form was mainly magnetically supported, and thus sub-Keplerian. Second, the 
substellar objects could transport some of their orbital angular momentum to the residual 
envelope by raising tides in it. In addition, since the disk contains more than one Jeans 
mass, more than one object could form from its fragmentation, with the ejection of all but 
one as a likely outcome of gravitational interactions. The final orbit of the remaining object 
could tighten by a large factor and tends to be highly eccentric (Papaloizou & Terquem 

2001) . Under extreme conditions, these interactions could in principle explain the smaller 
separations of Gliese 86B (~ 19 AU; Els et al. 2001) and HR 7672B (~ 14 AU; Liu et al. 

2002) . To send a brown dwarf even closer in, say to the ~ 3 AU distance of HD 168443c 
(Marcy et al. 2001), would be even more difficult, consistent with the well known dearth 
of brown dwarf companions within ~ 4 AU of solar-type FGK stars — the so-called "brown 
dwarf desert". For the same reason, our scenario may produce (massive) extrasolar giant 
planets on relatively wide orbits of tens to hundreds of AUs, but unlikely those within a few 
AUs of their host stars detected through radial-velocity surveys. These close-in planets are 
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presumably formed in rotationally supported disks, possibly through the conventional core- 
nucleation mechanism. Their orbital evolution may be affected, however, by the presence of 
substellar companions at larger distances. 

The substellar companions can also increase their orbital separations after formation, 
through gravitational interaction with the Keplerian disks they enclose. If formed in binary 
(or multiple) systems or dense clusters, they would likely be ejected through dynamic inter- 
actions, given their relatively large formation distances from the central star (Reipurth & 
Clarke 2001), although some may remain bound on wider orbits. Such interactions could in 
principle explain the the few thousand AU separations of the brown dwarf companions found 
in the 2 MASS survey and follow-up observations, including Gl 570D, Gl 417B, Gl 584C, Gl 
337C, and possibly HD 89744B and Gl 618. IB (Kirkpatrick et al. 2001; Wilson et al. 2001). 
Those ejected could account for, at least in part, the isolated brown dwarfs and possibly the 
free-floating planetary mass objects that may have been uncovered in a Orionis (Zapatero 
Osorio et al. 2000) and elsewhere. 

An esthetically appealing feature of our scenario is that the brown dwarf companions 
are produced in essentially the same manner as the isolated low-mass stars envisioned in 
the standard picture, through the magnetic diffusion-driven fragmentation of a strongly 
magnetized, self-gravitating medium. Two common products of the fragmentation in the 
stellar case are circumstellar disks and binaries. By analogy, we expect disks and binaries to 
result from the fragmentation in the substellar case as well. A surrounding disk could explain 
the infrared excess observed in the young brown dwarf in the GG Tau quadruple system 
(White et al. 1999), and brown dwarf pairs have been detected (Basri 2000). Dynamical 
ejection from the 10 2 AU-scale distance from a Sun-like low-mass star should leave intact the 
inner part of the disk (within a few AUs of the brown dwarf) or tightly bound pairs. The 
former is consistent with the inference of disks around young isolated brown dwarfs (Muench 
et al. 2001), and the latter with the dearth of isolated brown dwarf binaries with separations 
greater than ~ 10 AU (Reid et al. 2001). 

I thank G. Laughlin, C. McKee, F. Nakamura, R. Nishi, C. Terquem and the referee for 
helpful comments. 
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